Concatenated dynamical decoupling in a solid-state spin bath 
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Concatenated dynamical decoupling (CDD) pulse sequences hold much promise as a strategy to 
mitigate decoherence in quantum information processing. It is important to investigate the actual 
performance of these dynamical decoupling strategies in real systems that are promising qubit 
candidates. In this Rapid Communication, we compute the echo decay of concatenations of the 
Hahn echo sequence for a solid-state electronic spin qubit in a nuclear spin bath using a cluster 
expansion technique. We find that each level of concatenation reverses the effect of successive 
levels of intrabath fluctuations. On the one hand, this advances CDD as a versatile and realistic 
decoupling strategy. On the other hand, this invalidates, as overly optimistic, results of the simple 
pair approximation used previously to study restoration, through CDD, of coherence lost to a 
mesoscopic spin bath. 
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Successful large-scale quantum computation which re- 
lies exlusively on quantum error correcting codes neces- 
sitates extremely weak quantum decoherence (< 10 -4 ); 
this is difficult, if not impossible, to achieve in the 
strongly coupled solid-state environment. This stringent 
requirement arises from the fact that the quantum er- 
ror correction (QEC) overhead grows exponentially (e.g., 
the number of auxiliary qubits needed for encoding) with 
the magnitude of decoherence. At a local level, decoher- 
ence can be mitigated by logically encoding qubits in 
a subspace that is less sensitive to noise and by apply- 
ing dynamical decoupling (DD) controls. To meet QEC 
requirements, it is likely necessary to combine logical 
qubit encoding and DD in solid-state systemsi 1 - We focus 
our study on DD, which may employ deterministicf^^ 
random^ or hybrid deterministic-random^ sequence of 
qubit controls to partially reverse the effects of deco- 
herence. Deterministic sequences hold promise for slow, 
predictable bath systems such as the nuclear spin bath 
that we consider. In a universal DD sequence, the time- 
averaged Hamiltonian decouples the system from the 
bath, destroying the lowest-order system-bath coupling 
in a Magnus- expansion; 4 . The Hahn echo, with a sin- 
gle 7T pulse applied midway through the evolution, is a 
simple universal DD sequence for a dephasing qubit. Pe- 
riodic DD (PDD), repetitive applications of universal DD 
proven effective for nuclear magnetic resonance purposes 
for many years p simply trades an increase in the num- 
ber of applied pulses for an increase in their frequency. 
Higher orders of the Magnus expansion, however, limit 
the performance of PDD; it is possible to overcome this 
effect using concatenated DD (CDD), recursive applica- 
tions of a universal DD sequenced 

In this Rapid Communication, we consider an 
archetypical qubit, a localized electron spin in a semi- 
conductor at low temperature. The energy relaxation 
process of this spin can be made arbitrarily slow by low- 
ering the temperature (e.g., T\ ^> 1 ms for T ~ 100 mK) 
because phonons are exponentially suppressed below the 



Bloch-Gruneisen temperature (~ 10 K). In the limit of 
a strong applied magnetic field, slow fluctuations of the 
surrounding nuclear spins will induce electron spin de- 
phasing ("spectral diffusion"—) as the primary decoher- 
ence mechanism on a time scale of Ti ~ 1 fis (GaAs) — 
100 (is (Si:P). We show that this decoherence can be 
corrected, in principle, with arbitrary precision by us- 
ing concatenations of the Hahn echo. The true reason 
for this, interestingly, goes beyond consideration of or- 
ders in the Magnus expansion; each concatenation of the 
sequence reverses the effect of successive levels of intra- 
bath processes (specifically, CDD cancels orders of an 
intrabath perturbation). 

There have been general theoretical discussions 4 about 
using CDD for qubit preservation without quantitative 
applications to specific realistic models of qubits in meso- 
scopic baths. To the extent that specific mesoscopic ap- 
plications 1 ^ exist for the spin qubit example considered 
in this paper, they invariably involve highly simplified 
and uncontrolled (as well as untested) models, e.g., the 
nuclear spin pair approximation, which fails for the real- 
istic spin bath in the context of CDD. The pair approx- 
imation fails precisely because CDD, beyond canceling 
orders of a time (effectively Magnus) expansion, reverses 
the effect of successive levels of intrabath processes. As 
these lower levels of interactions are reversed (e.g., a sin- 
gle flip-flop of a pair of nuclei), larger clusters become 
perturbative equals to pairs and can dominate decoher- 
ence when they are greater in number. These larger clus- 
ters will diminish the effect of CDD (Refs. [13 give overly 
optimistic predictions). Nonetheless, we show that, con- 
trary to PDD, concatenation allows one to obtain the 
same fidelity while decreasing the frequency of applied 
pulses. Our work is, to the best of our knowledge, the 
first theoretical application of the CDD scheme for a re- 
alistic situation, establishing, as a matter of principle, 
the reasonable prospects for perpetual coherent control 
of a qubit. In proof of concept, we demonstrate this in 
the form of quantum memory retention; however, it is 
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known^ to be possible to effect nontrivial quantum evo- 
lution while performing DD. In fact, with proper logic 
encoding such that logical and DD operations commute, 
it may be possible^ to perform CDD and quantum logic 
simultaneously, making CDD a viable strategy to meet 
the error correction threshold for a fault-tolerant solid- 
state spin quantum computer. 

The continuous electron spin dephasing due to its cou- 
pling to the nuclear spin bath is a difficult process to ana- 
lyze theoretically, in spite of a very long history^ due to 
the quantum, non-Markovian nature of nuclear spin flip- 
flops. We previousl y 13 i 14 developed a quantum theory of 
nuclear-induced spectral diffusion using a formally exact 
cluster expansion technique and analyzed the Hahn echo 
decay of solid-state spin qubits, and later— adapted the 
technique to study the periodic Carr-Purcell-Meiboom- 
Gill 8 (CPMG) sequence. Concatenations of the Hahn 
echo sequence were analyzed in Refs. 10 for mesoscopic 
quantum-dot baths, which, however, employed a pair ap- 
proximation that is not accurate in light of the pertur- 
bative cancellations that CDD sequences induce. 

Specifically, we treat a donor electron spin in Si:P and 
a quantum dot electron spin in GaAs. Our electron spin 
qubit will precess in an applied magnetic field (~ 1 T) 
at its Zeeman frequency (~ 10 GHz) with some fluctu- 
ations due to hyperfmc (hf) interactions (~ 100 kHz) 
with a nuclear spin bath. The strong electron Zeeman 
energies compared with hf and nuclear Zeeman ener- 
gies (~ 10 MHz) suppresses direct hf interactions that 
would flip the electron spin. Most importantly, the nu- 
clei also interact with each other magnetically through 
their dipoles (~ 10 Hz) which, if unchecked, will lead to 
an irreversible loss of information. 

A general dephasing Hamiltonian may be written in 
the form H — ^ ± |±)7i±(±|, where fi± acts only upon 

the bath's Hilbert space. We can split H± into qubit- 
dependent and -independent parts, so that 7i± = ±7i q b + 
eHbb, where Ti. q b contains the qubit-bath interaction 
terms and Jibb contains the intrabath interaction terms. 
The e is a bookkeeping parameter that will be discussed 
later in terms of an intrabath perturbation. From the 
Fermi-contact hf interaction, H q b — \ '^Z n A n I nz , where 
A n is proportional to the probability of the electron being 
at the nth nuclear site. From the secular dipolar inter- 
action, 7~Lbb = y) n -t. m b nm I n +I m — — y] n yL m 2b nm I nz I mz 
where the first summation is restricted to pairs of like 
nuclei (so that Zeeman energy is preserved when they 
flip-flop) and b nm are dipolar coupling constants. Wc 
have included the most relevant interactions in the limit 
of a strong applied field as discussed in detail in Ref. [Til- 
Given an H = 5^ ± |±)H±(±|, the evolution opera- 
tor U for any sequence of tt pulses (acting perpendic- 
ularly to the applied field) may be written in the form 
U = ^ ± |±)C/ ± (±| (or U = E±lT)C> ± (±|) for an even 
(or odd) number of applied pulses. With the up-down 
free evolution operators denoted U^(r) = exp(i1-i±r), 
the Hahn echo sequence ofr— >7t^t^7t, with an 



extra ir pulse at the end to return the electron spin to its 
original orientation, has up-down evolution operators of 
Ut = U^{t)Uq{t). This trivially modified Hahn echo 
is the first level of the CDD series we consider. The 
Hahn echo sequence is best known as a means 1 ^ to elim- 
inate the effects of inhomogeneous broadening (reversing 
the effects of a distribution of electron Zeeman frequen- 
cies in an ensemble). However, the Hahn echo is also 
an elementary DD sequence because the time-averaged 
Hamiltonian, proportional to 7i + + Tt- oc Hbb, involves 
no qubit-bath interactions (in this respect, it does not 
yield a proper T2 time consistent with the free induction 
decay of a single qubit 16 ). As a DD sequence, it may be 
concatenated in hopes of successively improving qubit 
preservation times^ At the Zth level, our CDD pulse se- 
quence for I > is p ; := P;_i — > tt — > P;_x tt, and 
p := t 4 . With each concatenation, we do to the previ- 
ous sequence what the Hahn echo does to free evolution 
and in this way obtain improved DD. This sequence may 
be simplified by noting that two 7r pulses in sequence do 
nothing. Thus, 

(Pj-i-t-Pj-! , odd; 
I -> Pj-i , even I 

and the up-down evolution operators at level I have the 
recursive form of-£ 

= UT-M-v (2) 

In order to characterize the coherence decay, we con- 
sider the transverse component of the qubit 's expecta- 
tion value; normalized to a maximum of 1, the pulse se- 
quence echo ve is defined in this way such that ve — 
||([l>-] t L>+}|| = \\{W)\\ where W = [U~^U + and the 
(...) denotes an appropriately weighted average over the 
bath states. We assume that the bath is fully random 
and use equal weights in averaging over bath states^ 

We can solve this mesoscopic quantum problem using 
a cluster method i 13 i 14 Consider expanding W such that 
W = ^2n—oW[ n \ where contains contributions to 

W that involve n separate clusters of "operatively" in- 
teracting nuclei. To be specific, the set of nuclei involved 
in a term of must all be connected together via 

factors of bilinear interaction operators to form a single 
connected cluster. Clusters have spatial proximity when 
interactions are local. If it is possible to approximate 
(TyW) by including only clusters up to some small size 
that is much less than the number of nuclei in the bath, 
N, and if the bath is effectively uncorrelated initially 
(e.g., a random bath), then (W^) w (W^) n /n\. In this 
"cluster approximation," 

v E = \\(W}\\^exp(Re{(Wm}}). (3) 

This approximation can be tested by quantifying the con- 
tributions of extraneous terms due to overlapping clusters 
when distributing through (W^) n r^ Formally, (H'W} is 
the sum of all cluster contributions, (M^ 1 !) = 53c(^c)i 
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FIG. 1: Echo decay of an electron bound to a P donor in 
natural Si with an applied magnetic field along the [100] lat- 
tice direction for different DD pulse sequences as a function 
of the minimum time r between consecutive pulses. The dot- 
ted lines give corresponding results to the lowest order in the 
intrabath perturbation. Only the first 90% of the decay is 
shown to avoid displaying regions in which the cluster expan- 
sion fails to converge. Conveniently, coherence much better 
than 90% is necessary for QEC protocols and the long-time 
coherence decay is irrelevant. 

where the contribution from some cluster C may be com- 
puted using the recursive formula 

{W' e ) = {Wc)- £ HiW'J, (4) 

{Ci} disjoint, i 
Ci^0, dCC 

and (Wc) is simply the solution to (W) including only 
the nuclei in the set C. The cluster expansion results 
from expanding (M^ 1 !) in Eq. ([3]) to include clusters of 
successively increasing size. 

We make a comparison between our convergent clus- 
ter expansion results for CDD and earlier calculated 15 
CPMG pulse sequence echoes of both an electron bound 
to a P donor in natural Si (Fig. [1]) and a quantum dot 
electron in GaAs (Fig. [2]) plotted as a function of the 
minimum time t between consecutive pulses. The two- 
pulse CPMG sequence is the same as the CDD2 sequence 
(1 = 2 CDD), and longer even-pulsed CPMG sequences 
are simple repetitions of the CDD2 sequence. The no- 
table improvement of even-pulsed CPMG sequences over 
the Hahn echoi 5 - is directly due to its extra level of con- 
catenation. With an increasing number of CPMG pulses, 
the coherence as a function of t is diminished, but, as 
noted in Ref. [HI, the coherence as a function of the to- 
tal sequence time, t = Avt for a sequence with 2v pulses, 
tends to improve. With each level of concatenation in the 
CDD series, short-time coherence as a function of either 
t or the total sequence time t = 2 1+1 t improves. This 
demonstrates the restorative power of concatenation. 

We can understand the reason for CDD coherence 
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FIG. 2: Echo decay of a quantum dot (Fock-Darwin radius 
of I = 25 nm and quantum well thickness of zo = 8.5 nm) 
electron in GaAs with an applied magnetic field along the 
[110] lattice direction for different DD pulse sequences as a 
function of the minimum time r between consecutive pulses. 
The dotted lines give corresponding results to the lowest order 
in the intrabath perturbation. 

restoration as successive cancellations in perturbative or- 
ders. In a perturbation whose order increases with in- 
creasing cluster size, the lowest order of (W} — 1 is equiva- 
lent to the same order of (W^W) because terms of W with 
multiple clusters are automatically higher-order terms 
(being products of lower-order terms). This is formal 
equivalence that we use for convenience only (the en- 
tire purpose of the cluster expansion is to avoid applying 
such a perturbative expansion to (W) directly because it 
would fail to converge due to the large size of the bath). 

Using the unitarity property of the evolution operators 
in Wi = [U^Uf and defining A ; = U+ - U~ , 

Re {(#,)} = l-^A+A,), (5) 

so that Re{(Wi)} - 1 = -O (A?) = Re{(wf ] )} and 
Ve = 1 — O (A^J. Applying the recursive definitions for 
the (if 1 evolution operators [Eq. (|2J)] , 

A, = U+ Uf = [Uf_ v U+J = [Ur_ v Ai-i], (6) 

noting that Uf_ 1 commutes with itself. 

Consider a perturbation with a smallness parameter 
A in which Uf 1 = t + O (A) for all I > Iq for some ^o- 
Because the identity commutes with anything, Eq. ([6|) 
implies that A/ = O (A) x A/_i for all I > Iq. A pertur- 
bation in the time parameter r as well as a perturbation 
in the intrabath interactions e both have this property 
with Iq — and io = 1> respectively; CDD therefore 
cancels successive orders in r and e. The successive can- 
cellation in the time perturbation was previously noted 
in Refs. [lfj| . but not that in the intrabath perturbation. 
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FIG. 3: Contributions from clusters of different sizes for the 
concatenated echo decays in the (a) Si:P system of Fig. [1] 
and (b) GaAs system of Fig. [5] The minimum cluster size 
required to yield the appropriate lowest-order result in the 
cluster expansion increases with each concatenation, and the 
larger clusters tend to dominate the decoherence. 



Both of these perturbations are in correspondence with 
the cluster expansion because an increase in cluster size 
necessitates an increase in the minimum number of bi- 
linear interaction factors e, and a corresponding increase 
in t factors (by time-energy dimensional arguments). In 
fact, the cluster expansion converges due to the e pertur- 
bation (the time perturbation is additionally applicable 
only in special cases) as demonstrated in Figs. [1] and [H 
where the dotted lines give the results with (V^W) ap- 
proximated by the lowest order in e. This suggests that 
CDD provides relevant perturbative cancellations in the 
same regime where the cluster expansion converges. 

With each perturbative cancellation, clusters of in- 
creasing size must be included in the lowest-order clus- 



ter expansion calculation as all have the same pertur- 
bative order. Contributions from different cluster sizes 
are shown for the CDD series in Fig. [3] for the Si:P and 
GaAs systems. The larger clusters of the same perturba- 
tive order are seen to dominate the decoherence because 
they are greater in number (concatenation of level I is 
dominated by clusters of size I + 1 ) , invalidating the pair 
approximation used in Refs.fiol. This is unfortunate since 
the pair approximation is simple and physically transpar- 
ent. 

In summary, we establish the effectiveness of CDD 
in correcting spectral-diffusion-induced quantum deco- 
herence in the central spin problem. For our problem 
(and assuming ideal pulses), CDD effectively restores 
phase coherence when the interpulse delay is short com- 
pared with the typical intrabath interaction time scale. 
Furthermore, each level of concatenation involves larger 
clusters of bath spins, thus invalidating the simple pair 
approximation used earlier in the literature. In reality, 
the deviation from pulse ideality will determine the level 
of feasible quantum error correction; there is, however, 
promising recent work-i^ to address this issue. 

Our study of spin bath decoherence provides impor- 
tant tests for CDD strategies and takes these strategies 
beyond the level of abstract formalism 3 ' - and small-scale 
models^^ establishing that CDD could play a key role in 
correcting decoherence-induced errors in solid-state spin 
quantum computer architectures. An ability to prolong 
quantum coherence indefinitely using CDD techniques 
could lead to an effective solid-state spin quantum mem- 
ory as well as much reduced overhead in quantum error 
correction protocols. 
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